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Dynamical evolution plays a key role in shaping the current properties of star clus- 
ters and star cluster systems. A detailed understanding of the effects of evolutionary 
processes is essential to be able to disentangle the properties which result from dy- 
namical evolution from those imprinted at the time of cluster formation. In this 
review, we focus our attention on globular clusters and review the main physical 
ingredients driving their early and long-term evolution, describe the possible evo- 
lutionary routes and show how cluster structure and stellar content are affected by 
dynamical evolution. 
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1. Introduction 

Globular star clusters have long been considered the ideal astrophysical objects to 
explore many aspects of stellar dynamics and, in particular, to study the evolution 
of stellar systems as driven by the effects of two-body relaxation. 

Only recently, however, the actual complexity of globular cluster dynamics has 
emerged from the wealth of new observational and theoretical studies that have 
clearly shown the close interplay between stellar dynamics, stellar evolution, the 
clusters' stellar content and the dynamics and properties of the host galaxy. The 
notion of an 'ecology of star clusters', introduced by Heggie (1992), nicely illustrates 
the close interplay among the different elements of the astrophysics of star clusters. 

To use star clusters as tools to advance and guide our understanding of star for- 
mation in galaxies, it is essential to understand to what extent the current properties 
of star clusters and star cluster systems are the result of evolutionary processes. 
We also need to be able to discern the signatures of properties imprinted by for- 
mation processes from those determined by dynamical evolution. In this paper, we 
will review the main physical ingredients driving the early and long-term evolution 
of clusters and describe how cluster structure and stellar content is affected by 
dynamical evolution. 

In we discuss the early evolution of clusters, i.e., the processes that can 
affect their stellar content early in their life and lead to their rapid dissolution. In 

we review the long-term evolution of clusters as driven, primarily, by two-body 
relaxation, focusing our attention on the evolution of cluster structure, mass and 
stellar content. Conclusions are summarized in 521 

2. Early dynamical evolution 

Although much progress has been made in recent years (see, e.g., Clarke 2010; Lada 
2010), a full understanding of the formation and the initial structural properties of 
star clusters is still lacking. 
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Analytical and numerical studies of globular cluster evolution often adopt initial 
conditions which resemble the current structural properties of clusters (e.g., purely 
stellar, spherical and isotropic King models), rather than the possibly more complex 
structures suggested by theoretical and observational studies of very young star 
clusters (e.g., Elmegreen 2000; Bonnell et al. 2003; Allen et al. 2007). Figure 1 
shows a few snapshots of a hydrodynamical simulation modelling the formation 
of a small star cluster (Bonnell et al. 2003) and nicely illustrates the clumpy and 
irregular structure of a young star cluster still embedded in the primordial gas from 
which its stars are forming. 

In addition, the standard paradigm according to which clusters are 'simple stel- 
lar populations' composed of stars of the same age and chemical composition has 
recently been challenged by the increasing observational evidence of the presence of 
multiple stellar populations in globular clusters (e.g., Bedin et al. 2004; Gratton et 
al. 2004; Piotto et al. 2007; Carretta et al. 2008, 2009a, 6; D'Antona & Caloi 2008; 
Milone et al. 2008; see also the reviews of van Loon 2010; Kalirai & Richer 2010; 
and Bruzual 2010). This is a major paradigm shift and significantly increases the 
complexity of models following the formation and dynamical evolution of globular 
clusters (e.g., D'Ercole et al. 2008; Renzini 2008; and references therein). 

Although the exact structural and kinematic properties of newly formed clusters 
are still uncertain, a number of studies have explored the early evolution of clus- 
ters and shown their possible evolutionary routes. In the following subsections, we 
present an overview of the early dynamical evolution of clusters, based on analytical 
calculations and A-body simulations. 




Figure 1. Formation and early evolution of a star cluster from a hydrodynamical simulation 
following the hierarchical fragmentation of a turbulent molecular cloud (from Bonnell et 
al. 2003). 
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(a) Evolution of clumpy systems 



A number of observational and numerical studies suggest that star clusters might 
form with irregular clumpy structures and thus out of virial equilibrium. These 
studies mostly concern low-mass clusters and it is not clear whether such initial 
conditions are shared by more massive globular-cluster-like objects. Many clusters 
as young as <~ 50 — 100 Myr already exhibit smooth surface-brightness profiles 
(see, e.g., Elson et al. 1987; Larsen 2004) which are well fit by the 'Elson, Fall and 
Freeman' profile (hereafter EFF; Elson et al. 1987), 



where E is the cluster's surface brightness and r c its core radius. A number of studies 
have investigated the dynamics of clusters starting from clumpy initial conditions 
and explored if, and on what timescale, early dynamical evolution can drive them 
towards properties similar to those of older clusters with a smoother density profile. 

The evolution of clusters with different levels of dumpiness, clump spatial distri- 
butions and initial velocity distributions has been studied in general surveys (e.g., 
Aarseth et al. 1998; Goodwin 1998; Goodwin & Whitworth 2004; McMillan et al. 
2007) as well as in investigations aimed at modelling the evolution of specific clus- 
ters (e.g., Scally & Clarke 2002). Recently, Hurley & Bekki (2008) made further 
progress towards more realistic models and explored the stellar dynamical evolu- 
tion of an inhomogeneous, clumpy stellar system with initial conditions produced 
by an hydrodynamical simulation of a turbulent giant molecular cloud. 

The results of all these numerical studies show that, in most cases, an initially 
clumpy system undergoes a rapid phase of violent relaxation, which erases any 
initial substructure and leaves the cluster with a density profile which is very similar 
to the EFF profile. Similar simulations have also been carried out in the context of 
early galaxy evolution (sec, e.g., the first simulations by van Albada 1982; McGlynn 
1984; see Trenti et al. 2005, and references therein, for a more recent study) and, 
in fact, even before the EFF profile was used to fit the observed surface-brightness 
profiles of young clusters in the Large Magellanic Cloud, it was adopted by McGlynn 
(1984) in a study of early galaxy evolution to fit the density profiles emerging at 
the end of the violent-relaxation phase of initially clumpy systems. 

The timescale to reach final equilibrium without any significant remaining trace 
of the initial substructures depends on the details of the initial conditions but it 
is, in general, « 1 — 10 tdyn where £<jyn = /l&Gp is the cluster dynamical time 
(and p the cluster's average density). Fellhauer et al. (2009) recently showed that 
the clump-merging timescale may be shorter than the timescale of primordial-gas 
expulsion, thus enhancing the cluster's chances to survive the disruptive effects of 
gas expulsion (see §2 c below) 

Statistical tools to more rigorously measure cluster substructure have been used 
by Cartwright & Whitworth (2004), Schmeja & Klessen (2006) and Kumar & 
Schmeja (2007) and provided a quantitative measure of the evolving level of dumpi- 
ness in observed young clusters and star- forming regions, as well as in the numerical 
results of hydrodynamical simulations of star cluster formation. 
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(b) Clumpy systems and initial mass segregation 



Although the evolution of clumpy stellar systems is rapid and occurs on timescales 
on the order of a few dynamical timescales, it can have important implications for 
the generation of mass segregation early in a cluster's evolution. Mass segregation, 
the tendency of more massive stars to sink towards the cluster's central regions, is 
a consequence of energy equipartition driven by two-body relaxation and one of the 
processes normally associated with the long-term evolution of clusters. However, a 
number of observational studies (e.g., Hillenbrand 1997; Hillenbrand & Hartmann 
1998; Fischer et al. 1998; de Grijs et al. 2002; Sirianni et al. 2002; Gouliermis et 
al. 2004; Stolte et al. 2006; Sabbi et al. 2008; Allison et al. 2009a; Gouliermis et 
al. 2009; see also Ascenso et al. 2009) have found evidence of mass segregation in 
clusters with ages shorter than the time needed to produce the observed segregation 
by two-body relaxation (see also de Grijs 2010). 

It has been suggested in a number of theoretical studies that the observed seg- 
regation in young clusters would be primordial and imprinted by the star-formation 
process (Bonnell et al. 1997, 2001; Bonncll & Davies 1998; Klessen 2001; Bonncll & 
Bate 2006). Specifically, the mechanism responsible for primordial mass segregation 
would be the higher accretion rate for stars in the central regions of young clus- 
ters. The actual efficiency of this mechanism is still a matter of debate (Krumholz 
et al. 2005; Krumholz & Bonnell 2007), but if the individual clumps are indeed 
mass segregated, it has been shown by McMillan et al. (2007) that such primordial 
mass segregation would not be erased in the violent-relaxation phase during which 
clumps merge. The final system would preserve the mass segregation of the original 
clumps (see also Fellhauer et al. 2009; Moeckel & Bonnell 2009). 

In the same study, McMillan et al. also presented an alternative scenario for a 
dynamical origin of early mass segregation in young clusters. Even if the clumps are 
not initially segregated, if their internal segregation timescale is shorter than the 
time needed for the clumps to merge, they will segregate through standard two-body 
relaxation and preserve this segregation after they have merged. The multiscale dy- 
namical evolution of clumpy systems is, in this case, responsible for rapidly leading 
to mass segregation in young clusters without invoking any mechanism associated 
with the star-formation process. 

In a recent study, Allison et al. (20096) showed that mass segregation can be 
rapidly produced dynamically also in the high-density core formed during the col- 
lapse of cold clumpy clusters. Bastian et al. (2008) found observational evidence of 
a strong expansion in the first 20 Myr of evolution of six young M51 clusters and 
pointed out that this expansion could also lead to a rapid variation in the cluster 
relaxation time. As pointed out by Bastian et al, using the current relaxation time 
might lead to an underestimate of the possible role played by two-body relaxation in 
generating mass segregation in the early phases of a cluster's dynamical evolution. 

Regardless of the mechanism producing mass segregation, the presence of segre- 
gation very early on in a cluster's life can have a significant impact on its dynamical 
evolution. We will discuss some of the implications of initial mass segregation in 
the next subsection. 
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(c) Early mass loss: primordial gas expulsion and stellar evolution 

Expulsion of the primordial residual gas in which a young cluster is still embed- 
ded in the very early stages of its life, along with mass loss due to stellar evolution, 
can have a significant impact on a cluster's structure and survival chances. 

Early analytical calculations by Hills (1980) based on the virial theorem showed 
that clusters in virial equilibrium losing impulsively more than half of their mass 
would rapidly expand and dissolve in response to this mass loss. Subsequent studies 
by Boily & Kroupa (2003a ; &) refined Hills' results by calculating the fraction of stars 
remaining bound in a cluster with a given velocity distribution function after the 
impulsive removal of a given amount of gas. Boily & Kroupa's study showed that, 
for clusters with a high-binding-energy massive core, up to 70% of the cluster mass 
can be removed impulsively without leading to complete cluster dissolution (see 
also the simulations in Kroupa et al. 2001; and references therein). 

Goodwin & Bastian (2006) explored, by means of A^-body simulations, the evo- 
lution of the structural properties of clusters in the stages following gas expulsion 
and showed that the observed luminosity profiles of a few young massive clusters 
differ from EFF and King models in their outer regions, consistent with their sim- 
ulations. These could, therefore, be examples of clusters currently expanding in 
response to gas expulsion. 

Baumgardt & Kroupa (2007) carried out a survey of iV-body simulations to ex- 
plore the dependence of the response to residual gas expulsion on the star-formation 
efficiency, the expulsion timescale and the strength of the host galaxy's tidal field. 
Slow gas removal and a weak tidal field can increase the amount of mass a cluster 
can lose without undergoing complete disruption. However, in most cases a cluster 
undergoes significant expansion and needs to be initially much more compact than 
observed today to survive this phase. As shown in a number of numerical studies 
(e.g., Portegies Zwart et al. 2002, 2004; Gfirkan et al. 2004; see also McMillan 2008, 
and references therein, for a review), very high initial central densities might lead 
to a rapid core collapse and segregation of massive stars, and trigger a runaway 
merger of massive stars, leading to the formation of an intermediate-mass black 
hole (IMBH; but see Glebbeek et al. 2009 for a recent study of the stellar evolution 
of runaway merger products, showing how mass loss due to stellar winds might 
significantly affect this process and prevent the formation of an IMBH). 

Early mass loss due to stellar evolution (for example, expulsion of Type II su- 
pernova ejecta) can also have a significant impact on early cluster evolution. Semi- 
analytic calculations (e.g., Applegate 1986; Chernoff & Shapiro 1987) followed by 
Fokker-Planck (e.g., Chernoff & Weinberg 1990) and A^-body simulations (e.g., 
Fukushigc & Heggie 1995; Portegies Zwart et al. 1998) have shown that mass loss 
due to stellar evolution triggers an expansion, which can lead to rapid cluster dis- 
solution (see also Kalirai & Richer 2010). 

The extended surveys of Chernoff & Shapiro (1987) and Chernoff & Weinberg 
(1990) showed that mass loss due to stellar evolution can cause the rapid dissolution 
of clusters with a low initial concentration and/or a flatter stellar initial mass 
function (IMF). 

Fukushige & Heggie (1995) used A^-body simulations to explore the early evolu- 
tion of clusters dissolving due to mass loss associated with stellar evolution. Their 
results are in general, qualitative agreement with those of the Fokker-Planck mod- 
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els of Chcrnoff & Weinberg (see also Takahashi & Portegies Zwart 2000 for an 
extensive comparison between TV-body simulations and Fokker-Planck models). 
Fukushige & Hcggic also explored the mechanism behind this rapid cluster disso- 
lution and showed it to be the result of a loss of equilibrium as the cluster expands 
and reach structural properties for which no virial equilibrium is possible. 

The expansion triggered by early mass loss due to stellar evolution is stronger 
if the cluster is initially mass segregated. For a given amount of mass loss, pref- 
erentially removing mass from the central regions (where massive stars tend to 
be located in mass-segregated clusters) increases the heating and strengthens the 
subsequent expansion. TV-body simulations by Vesperini et al. (2009) show that as 
the degree of initial mass segregation increases, so does the strength of the initial 
cluster expansion. As a result, clusters differing only in the degree of initial mass 
segregation can have very different lifetimes. Mackey et al. (2007, 2008) showed that 
the stronger early expansion of mass-segregated clusters, along with the subsequent 
heating from a population of stellar black holes, can explain the radius-age trend 
observed for massive clusters in the Magellanic Clouds. 

A strong early loss of stars plays a key role in the evolution of multiple stellar 
populations in clusters. Although different origins have been suggested for the gas 
from which second-generation stars might have formed (e.g., Decressin et al. 2001; 
Ventura et al. 2001; see also Renzini 2008; and references therein), most models 
require the first-generation population to have an initial mass which is at least ten 
times higher than its current mass. As shown by D'Ercole et al. (2008), early cluster 
expansion and mass loss can be responsible for the escape of such a large fraction 
of first-generation stars. 

Finally, although this review focuses on the internal dynamics of clusters rather 
than on cluster systems (see Harris 2010; Larsen 2010 in this volume), it is important 
to point out that evolution and disruption of clusters due to gas expulsion and 
stellar evolution may play an important role in the evolution of the mass function 
of globular cluster systems (hereafter GCMF). 

Young cluster systems have a power-law GCMF (e.g., Zhang & Fall 1999) while 
older clusters follow a bell-shaped distribution (e.g., Ashman & Zepf 1998; Harris 
2001; Brodie & Strader 2006). While evaporation due to two-body relaxation can 
transform a power law into a bell-shaped GCMF, most GCMF evolution models 
show that when GCMF evolution is driven only by two-body relaxation, a galacto- 
centric dependence of the GCMF turnover not found in observational data is also 
produced (e.g., Vesperini 1997, 1998, 2000, 2001; Baumgardt 1998; Fall & Zhang 
2001; Vesperini et al. 2003; McLaughlin & Fall 2008; Baumgardt et al. 2008 a; see 
Zepf 2008 for a recent review). 

It has been shown that early cluster dissolution due to gas expulsion and mass 
loss from stellar evolution can preferentially destroy low-mass clusters and signifi- 
cantly flatten the low-mass end of the initial power-law GCMF, without introducing 
any strong dependence on galactocentric distance. Two-body relaxation and tidal 
shocks, while still playing a key role and leading to disruption of a significant 
number of clusters, will act on the GCMF already flattened by early evolutionary 
processes without giving rise to any radial trend of GCMF properties, inconsistent 
with observations (Vesperini & Zepf 2003; Baumgardt et al. 2008 a; Parmentier et 
al. 2008). 
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3. Long-term evolution of clusters 

For most old globular clusters, the core and half- mass two-body relaxation timescales 
are shorter than the cluster ages. Two-body relaxation plays a major role in driving 
the long-term evolution of clusters and shaping their current structural properties, 
as well as their stellar content. 

In this section, we focus our attention on the effects of two-body relaxation. We 
point out, however, that the evolution of clusters on eccentric orbits and clusters in 
spiral galaxies can be further affected by the tidal shocks during passages through 
the disc or near the host galaxy's central bulge. Although in most cases two-body 
relaxation is the dominant evolutionary process, depending on the properties of the 
host galaxy and the cluster's orbital parameters, tidal shocks can significantly speed 
up cluster dissolution and either accelerate or slow down the evolution towards core 
collapse (e.g., Chernoff et al. 1986; Spitzer 1987; Weinberg 1994; Gnedin et al. 1999). 
Encounters with giant molecular clouds and passages through spiral arms can also 
affect cluster evolution and mass loss (Gieles et al. 2006, 2007). 

In the following subsections, we focus our attention on mass loss due to two- 
body relaxation, the implications of mass loss for the evolution of the cluster stellar 
mass function and the cluster's structural evolution. 

(a) Mass loss and cluster dissolution 

As a star exchanges energy with other single and binary stars in a cluster due to 
close and distant encounters, it can reach an energy in excess of the escape energy 
and escape from the cluster. 

Mass loss due to two-body relaxation plays a key role in a broad range of 
issues related to the evolution of globular clusters and globular cluster systems. 
Specifically, for clusters surviving the early evolutionary processes described in fj2j 
mass loss due to two-body relaxation is the main process determining their lifetimes. 
It affects the clusters' stellar content and stellar mass function and determines the 
fraction of stars in the host galaxy's field population contributed by clusters. 

Although mass loss due to two-body relaxation has been the subject of a large 
number of studies since the very early investigations of star cluster dynamics (see, 
e.g., Heggie 2001 for a review, and references therein), only recently a number of 
investigations have shed light on some fundamental aspects of this process and its 
dependence on cluster structural parameters and the strength of the host galaxy's 
tidal field. 

The external tidal field of the host galaxy plays a key role in the process of 
mass loss. While isolated clusters undergo mass loss due to the combined effects of 
close and distant encounters (e.g., Giersz & Heggie 1994; Heggie 2001; Baumgardt 
et al. 2002; Heggie & Hut 2003), their dissolution timescale is extremely long. The 
N-body simulations of isolated star clusters of Baumgardt et al. (2002) show that 
it takes about 10 3 i r h(0) for a cluster to lose about half of its initial mass, where 
i r h(0) is the cluster's initial half-mass relaxation time. 

For clusters evolving in the external tidal field of the host galaxy, on the other 
hand, the mass-loss rate is much faster and the dissolution time significantly shorter: 
by lowering the escape speed and truncating the cluster sizes, the external tidal 
field significantly enhances the mass-loss rate. The process of mass loss due to two- 
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body relaxation and the dissolution time associated with this process depend on 
the strength of the external tidal field. Fokker-Planck (e.g., Chernoff & Weinberg 
1990; Takahashi & Portegies Zwart 2000) and A-body simulations (e.g., Vesperini 
& Heggie 1997; Aarseth & Heggie 1998) have shown that the mass-loss rate and 
dissolution time, T d , of a cluster containing an initial number of stars A;, evolving on 
a circular orbit at a galactocentric distance R g in a host galaxy with circular velocity 
v c is proportional to i og ( N .^ v ■ For given values of Aj, w c and R g , the mass-loss rate 
depends weakly on the cluster's internal structural properties (e.g., concentration 
and half-mass radius). An extensive survey of TV-body simulations was recently 
carried out by Gieles & Baumgardt (2008). Their results clearly show that, for 
a cluster with a given Ai, the strength of the tidal field is the dominant factor 
determining a cluster's lifetime. 

If the presence of the tidal field is properly modelled (as opposed to being 
treated in a simplified way by introducing, for example, a spatial or an energy 
cutoff), stars can escape the cluster only through one of the Lagrangian points 
of the galaxy-cluster system. Fukushige & Heggie (2000) and Baumgardt (2001) 
showed that the time needed for a star with an energy greater than the escape 
energy to flow through one of the Lagrangian points and actually escape from the 
cluster is not negligible and can, depending on the star's orbital parameters, be 
even longer than a Hubble time. Fukushige & Heggie (2000), Baumgardt (2001) 
and Baumgardt & Makino (2003) showed that the escape timescale modifies the 
scaling of the cluster dissolution time with cluster mass, and that Td is actually 
proportional to (Ni/log(Ni)) ' 75 (see also Lamers et al. 2005 for analytical fits to 
the mass-loss rates resulting from the A-body simulations of Baumgardt & Makino 
2003). 

A number of studies based on A-body simulations have followed the motion of 
escaping stars beyond the cluster's tidal radius and explored the formation, struc- 
ture and stellar content of the elongated tidal tails that these stars populate and 
the relationship between their orientation and the cluster's orbit (e.g., Combes et 
al. 1999; Johnston et al. 1999; Dehncn et al. 2004; Fellhauer et al. 2007; Mon- 
tuori et al. 2007; Ernst et al. 2009; Odenkirchcn et al. 2009). Observational studies 
have detected the presence of tidal tails in several clusters (e.g., Leon et al. 2000; 
Odenkirchcn et al. 2003; Grillmair & Dionatos 2006) and also revealed the presence 
of density clumps in the tidal tails (e.g., Leon et al. 2000). Similar clumps were 
found by Capuzzo Dolcetta et al. (2005) in A-body simulations of clusters on ec- 
centric orbits. More recently, Kupper et al. (2008) and Just et al. (2009) carried out 
detailed studies of the orbits of escaping stars and showed that overdensities along 
the tails are expected even for clusters on circular orbits. The observed clumps can 
be explained in terms of the epicyclic motion of escaping stars and a slowdown in 
some portions of their orbits as they move away from the cluster. 

As we will discuss in the next subsection, many theoretical studies have shown 
that mass loss due to two-body relaxation leads to the preferential escape of low- 
mass stars. In agreement with these results, Koch et al. (2004) — in an observational 
study of the Pal 5 main-sequence stellar luminosity function — found a significant 
overabundance of low-mass stars in this cluster's extended tidal tails. 
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(b) Mass loss and evolution of the stellar mass function 

As a result of the tendency towards energy equipartition, massive stars tend to 
sink to the cluster's central regions, while low-mass stars populate the outer parts 
and dominate the population of escaping stars. 

A consequence of the preferential loss of low-mass stars is that the stellar mass 
function flattens as a cluster proceeds in its dynamical evolution and loses mass. 
A number of iV-body simulations have revealed a close correlation between the 
fraction of the initial cluster mass lost and the slope of the low-mass end of the 
stellar mass function (Vesperini & Heggie 1997; Baumgardt & Makino 2003; Trenti 
et al. 2009). Assuming that clusters are characterized by a universal stellar IMF, 
the current slope of the mass function is a good indicator of the extent to which a 
cluster has been affected by dynamical evolution. 

The relationship between a cluster's dynamical evolution, mass loss and its stel- 
lar mass function is one of the manifestations of the close interplay between different 
aspects of the astrophysics of star clusters. Observational investigations of the stel- 
lar mass function of a number of Galactic globular clusters have found a range of 
different slopes of the low-mass side of the present-day mass function, which might 
indeed be the result of differences in the clusters' dynamical histories (cf. Piotto & 
Zoccali 1999). A recent observational study by De Marchi et al. (2007) resulted in 
a puzzling trend between cluster structure and slope of the stellar mass function: 
clusters with a flatter present-day mass function tend to have low-concentration 
density profiles. The origin of this trend is not clear (see, e.g., Baumgardt et al. 
2008 & for a study suggesting that primordial mass segregation might be required 
to explain the extreme mass function flattening of some low-concentration clusters) 
and will require further investigation of the evolution of cluster structural properties 
as they evolve towards complete dissolution. 

A consequence of the preferential loss of low-mass stars is that the evolution 
of a cluster's mass-to-light ratio, M/L, differs from that driven solely by stellar 
evolution. Specifically, it has been shown by Baumgardt & Makino (2003; see also 
Lamers et al. 2006; Anders et al. 2009 for recent studies of the photometric evolution 
of dissolving clusters) that during most of a cluster's life, M/L is lower than what 
the cluster would have without any loss of stars. Only at the very end of its life, 
when the cluster approaches complete dissolution, its stellar population becomes 
dominated by dark remnants and M/L rises to values greater than those expected 
if there were no mass loss. Kruijssen (2008) and Kruijssen & Mieske (2009) further 
investigated the effect of the preferential loss of low-mass stars on M/L and showed 
that the effects of mass loss can account for the discrepancy between the observed 
M/L values for a sample of Galactic globular clusters and those expected on the 
basis of including only the effects of stellar evolution. 

Tidal shocks due to passages near the host galaxys's bulge for clusters on ec- 
centric orbits, and through the Galactic disc, will inject energy into the cluster and 
further speed up the process of mass loss due to two-body relaxation (e.g., Gncdin 
& Ostriker 1997; Gnedin et al. 1999). However, mass loss due to tidal shocks is in- 
dependent of stellar mass and only combined with mass segregation can mass loss 
due to tidal shocks contribute to the flattening of the stellar mass function (e.g., 
Vesperini & Heggie 1997). 
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Figure 2. Time evolution of the 0.1, 1, 10 and 50% Lagrangian radii and tidal radius from 
an iV-body simulation following the dynamical evolution of a star cluster orbiting at 4 kpc 
from the Galactic centre and with an initial mass of M = 1.49 X 10 5 M Q (from Aarseth 
& Heggie 1998). 

(c) Structural evolution, core collapse and post- core collapse 

Figure 2 shows the results of an TV-body simulation (from Aarseth & Heggie 
1998) following the time evolution of the 0.1, 1, 10 and 50% Lagrangian radii and 
tidal radius of a star cluster orbiting at 4 kpc from the Galactic centre and with 
an initial mass of Mj ~ 1.5 x 10 5 M©. The initial expansion of the cluster is a 
consequence of mass loss due to stellar evolution. As discussed in ^ clusters with 
low initial concentrations or a large fraction of massive stars, or with strong initial 
mass segregation, can undergo significant early expansion and rapidly dissolve as a 
consequence of mass loss due to stellar evolution. 

For clusters surviving early dissolution, figure 2 illustrates the typical structural 
evolution towards core collapse driven by two-body relaxation. First studied by 
Henon (1961) using a Monte Carlo integration of the Fokker-Planck equation, this 
process, also known as gravothermal catastrophe, has subsequently been explored 
in a large number of investigations (Antonov 1962 and Lynden-Bell & Wood 1968 
are among the first studies that explored the physics of this process; see Heggie & 
Hut 2003 for a review, and references therein). Without intervention by an energy 
source balancing the loss of energy from the core, this process would lead to a 
smaller and smaller core and a diverging central density in a finite time. 

Binary stars, either primordial or dynamically formed during close encounters 
between single stars, can provide the energy needed to halt core collapse and sup- 
port the core in the post-core-collapse (PCC) phase (see, e.g., Heggie & Hut 2003; 
and references therein, also for a review of the gravothermal oscillations that can 
characterize the PCC phase of clusters). In particular, as shown by Heggie (1975; 
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see also Heggie & Hut 1993, 2003 for reviews and Goodwin 2010), binary stars with 
binding energy, eb, moving in a stellar system of single stars with mass m and one- 
dimensional velocity dispersion a such that \e-t,\/ma 2 > 1 will, on average, increase 
their binding energy and release the energy lost to the star cluster. 

The depth of core collapse and the concentration of a cluster in the PCC phase, 
as measured, for example, by the ratio of the core to half-mass radius, r c /rh, depends 
on the energy source. Clusters supported by dynamically formed binaries will, in 
general, be extremely concentrated, r c /rh ~ 10~ 3 (e.g., Heggie & Hut 2003). 

The concentration of a cluster supported by primordial binaries depends on the 
abundance of binaries and, to a smaller extent, on the distribution of binary binding 
energy. Primordial binaries can support a cluster in the PCC phase with values of 
r c /rh as large as ~ 0.05 — 0.08 (see the analytical studies of Goodman & Hut 
1989; Vesperini & Chernoff 1994; and the Fokker-Planck and A^-body simulations 
of McMillan et al. 1990, 1991; Gao et al. 1991; Heggie & Aarseth 1992; Gicrsz & 
Spurzcm 2000, 2003; Heggie et al. 2006; Fregeau & Rasio 2007; Portegics Zwart et 
al. 2007; Trcnti et al. 2007a). 

Observational studies have shown that about 20% of Galactic clusters have 
cuspy surface-brightness profiles and identified these objects as clusters in the PCC 
phase (Djorgovski & King 1986; Chernoff & Djorgoski 1989; Djorgovski & Meylan 
1994). As predicted by Chernoff & Shapiro (1987), PCC clusters are preferentially 
located near the Galactic centre (Chernoff & Djorgovski 1989). The smaller cluster 
sizes set by the stronger tidal field at small galactocentric distances lead to shorter 
relaxation times and more rapid evolution towards core collapse. 

While many clusters are still in the pre-core-collapse phase, a few clusters have 
central and half-mass relaxation times which are short enough that they should 
have already undergone core collapse (Chernoff & Djorgovski 1989; Trcnti 2006). 
These clusters, however, have normal surface-brightness profiles, with a resolved 
core and values of r c /rh greater than can be supported even by a large fraction of 
primordial binaries (and it is important to note that recent observational studies 
suggest that the fraction of primordial binaries in several clusters might only be 
w 1 - 2%; Davis et al. 2008 a) . 

Alternative energy sources have been suggested. Trenti (2006) suggested that 
an IMBH at the centre of these clusters might support the large values of r c /rh 
(sec Baumgardt et al. 2005; Heggie et al. 2007; Trenti et al. 20076; Gill et al. 2008; 
Pasquato et al. 2009 for A^-body simulations of clusters containing an IMBH). 
Heyl (2008), Heyl & Pcnrice (2009) and Fregeau et al. (2009), following recent 
observational studies suggesting white dwarfs might receive a kick velocity at the 
time of their formation (Davis et al. 20086; see also Kalirai & Richer 2010), have 
explored the effects of these kicks on cluster dynamics and showed that they can 
indeed represent a significant energy source capable of delaying core collapse and 
supporting high values of r c /rh. 

Hurley (2007) and Trenti et al. (2009) analysed the results of a number of N- 
body simulations, following observational procedures, and showed that adopting the 
observational definitions of clusters' structural parameters reduces the differences 
between observational and theoretical values of parameters such as the ratio of the 
core to half-mass radius. As numerical models become more realistic, adopting data- 
analysis procedures consistent with those used in observational studies is important 
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and will help to clarify for which clusters an additional energy source is indeed 
required to explain the current cluster structure. 

The results of detailed iV-body models for the dynamical evolution of M4 and 
NGC 6397 (Heggie & Giersz 2008, 2009 a, b) show that there are additional factors 
that can complicate the observational identification of a cluster's dynamical phase. 
M4 and NGC 6397 are two Galactic clusters with expected similar dynamical history 
which, however, have different surface-brightness profiles: M4 has a normal King 
profile (although of high concentration) while NGC 6397 exhibits a cuspy profile and 
is usually classified as a PCC cluster. According to the TV-body models of Heggie & 
Giersz (2008, 2009a, 6), both M4 and NGC 6397 are actually in the PCC phase and 
the differences in the surface-brightness profiles are due to oscillations in the core 
size which, in turn, cause the surface-brightness profile to oscillate between cuspy 
and high-concentration King profiles. Depending on the phase of this oscillatory 
behaviour, a cluster will be classified either as a normal, high-concentration King 
model or a PCC cluster. 

The results of all these studies suggest that it is likely that not all clusters 
in the PCC phase have been classified as PCC clusters. Further exploration of 
the range of possible observational properties of PCC clusters and of additional 
observational indicators of the PCC phase is necessary to ensure that the dynamical 
state of normal and PCC clusters is correctly identified and that a proper study of 
correlations between a cluster's dynamical state and other internal (e.g., abundance 
of exotic objects such as pulsars, blue stragglers, etc.) and external properties (e.g., 
position in the Galaxy) can be carried out. 

4. Conclusions 

The study of star cluster dynamics has a long history and significant progress 
has been made in understanding the physics of the main ingredients driving the 
dynamical evolution of star clusters. Many fundamental problems, however, are 
still open and new observations continuously present us with new challenges and 
questions. 

Understanding the role played by early evolutionary processes in the dissolution 
of star clusters and in shaping the properties of those which survive, how to identify 
a cluster's dynamical state, what energy sources support clusters in the PCC phase 
and what the role is of evolutionary processes in determining the current properties 
of globular cluster systems are just a few of the issues that are still awaiting firm 
solutions. 

As the close interplay between dynamical evolution, the evolution of a cluster's 
stellar content and the role played by the host galaxy's environment has become 
clear, new and increasingly complex questions crossing the boundaries of different 
fields have arisen. Understanding the links between the formation and abundance 
of exotic objects (such as blue stragglers, X-ray sources and IMBHs) and a clus- 
ter's dynamical history, the relationship between a cluster's stellar mass function, 
its structure and its past dynamical evolution, the role played by the host galaxy's 
tidal field and the effect of its time variation during galaxy assembly in shaping 
the current properties of individual clusters and the global properties of cluster 
populations are some examples of unsolved problems that require increasingly so- 
phisticated models. The observational evidence of the presence of multiple stellar 
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populations in globular clusters is the latest major challenge and further illustrates 
the complexity of the formation and dynamical evolution of star clusters. 

As clusters evolve, their structural properties and stellar content are modified by 
evolutionary processes. Clusters in different galaxies and at different galactocentric 
distances have different dynamical histories. A better understanding of the effects 
of dynamical evolution is an essential step in our attempts to shed light on the 
relationship between current star cluster properties and those which were imprinted 
by star-formation processes. 

I acknowledge partial research support from NASA grant NNX08AHf5G. 

References 

Aarseth, S. J., Lin, D. N. C. & Papaloizou, J. C. B. f988 On the collapse and violent 

relaxation of protoglobular clusters. Astrophys. J., 324, 288-310. 
Aarseth, S. J. & Heggie, D. C. 1998 Basic iV-body modelling of the evolution of globular 

clusters. I. Time scaling. Mon. Not. R. Astron. Soc, 297, 794-806. 
Allen, L., Megeath, S. T., Gutermuth, R., Myers, P. C, Wolk, S., Adams, F. C, Muzerolle, 

J., Young, E. & Pipher, J. L. 2007 The structure and evolution of young stellar clusters. 

In Protostars and planets V(eds B. Reipurth, D. Jewitt & K. Keil), pp. 361-376, Tucson: 

University of Arizona Press. 
Allison, R. J., Goodwin S. P., Parker R. J., Portegies Zwart S. F., de Grijs R. & Kouwen- 

hoven M. B. N. 2009a Using the minimum spanning tree to trace mass segregation. 

Mon. Not. R. Astron. Soc, 395, 1449-1454. 
Allison, R. J., Goodwin, S. P., Parker, R. J., de Grijs, R., Portegies Zwart, S. F. & 

Kouwenhoven, M. B. N. 2009b Dynamical mass segregation on a very short timescale. 

Astrophys. J., 700, L99-L103. 
Anders, P., Lamers, H. J. G. L. M. & Baumgardt, H. 2009 The photometric evolution of 

dissolving star clusters. II. Realistic models. Colours and M/L ratios. Astron. Astro- 
phys., 502, 817-832. 

Antonov, V. A. 1962 Most probable phase distribution in spherical star systems and 
conditions for its existence. Vest. Leningrad Univ., 7, 135. Transl. in Dynamics of Star 
Clusters, eds J. Goodman & P. Hut, pp. 525-540, Dordrecht: Reidel. 

Applegate, J. H. 1986 Dynamical effects of stellar evolution in globular clusters. Astrophys. 
J., 301, 132-144. 

Ascenso, J., Alves, J. & Lago, M. T. V. T. 2009 No evidence of mass segregation in massive 

young clusters. Astron. Astrophys., 495, 147-155. 
Ashman, K. M. & Zepf, S. E. 1998 Globular Cluster Systems. Cambridge: Cambridge 

University Press. 

Bastian, N., Gieles, M., Goodwin, S. P., Trancho, G., Smith, L. J., Konstantopoulos, I. 
& Efremov, Yu. 2008, The early expansion of cluster cores. Mon. Not. R. Astron. Soc., 
389, 223-230. 

Baumgardt, H. 1998 The initial distribution and evolution of globular cluster systems. 

Astron. Astrophys., 340, 402-414. 
Baumgardt, H. 2001 Scaling of iV-body calculations. Mon. Not. R. Astron. Soc, 325, 

1323-1331. 

Baumgardt, H., Hut, P. & Heggie, D. C. 2002 Long-term evolution of isolated iV-body 

systems. Mon. Not. R. Astron. Soc, 336, 1069-1081. 
Baumgardt, H. & Makino, J. 2003 Dynamical evolution of star clusters in tidal fields. 

Mon. Not. R. Astron. Soc, 340, 227-246. 



Article to appear in the Philosopical Transactions of the Royal Society A 



14 



E. Vesperini 



Baumgardt, H., Makino, J. & Hut P. 2005 Which globular clusters contain intermediate- 
mass black holes? Astrophys. J., 620, 238-243. 

Baumgardt, H. & Kroupa, P. 2007 A comprehensive set of simulations studying the in- 
fluence of gas expulsion on star cluster evolution. Mon. Not. R. Astron. Soc, 380, 
1589-1598. 

Baumgardt, H., Kroupa, P. & Parmentier, G. 2008a The influence of residual gas expulsion 
on the evolution of the Galactic globular cluster system and the origin of the Population 
II halo. Mon. Not. R. Astron. Soc, 384, 1231-1241. 

Baumgardt, H., De Marchi, G. & Kroupa, P. 20086 Evidence for primordial mass segre- 
gation in globular clusters. Astrophys. J., 685, 247-253. 

Bedin, L. R., Piotto, G., Anderson, J., Cassisi, S., King, I. R., Momany, Y. & Carraro, G. 
2004 u) Centauri: the population puzzle goes deeper. Astrophys. J., 605, L125-L128. 

Boily, C. M. & Kroupa, P. 2003a The impact of mass loss on star cluster formation. I. 
Analytical results. Mon. Not. R. Astron. Soc, 338, 665-672. 

Boily, C. M. & Kroupa, P. 20036 The impact of mass loss on star cluster formation. II. 
Numerical JV-body integration and further applications. Mon. Not. R. Astron. Soc, 
338, 673-686. 

Bonnell, I. A., Bate, M. R., Clarke, C. J. & Pringle, J.E. 1997 Accretion and the stellar 

mass spectrum in small clusters. Mon. Not. R. Astron. Soc, 285, 201-208. 
Bonnell, I. A. & Davies, M. B. 1998 Mass segregation in young stellar clusters. Mon. Not. 

R. Astron. Soc, 295, 691-698. 
Bonnell, I. A., Clarke, C. J., Bate, M. R. & Pringle, J. E. 2001 Accretion in stellar clusters 

and the initial mass function. Mon. Not. R. Astron. Soc, 324, 573-579. 
Bonnell, I. A., Bate, M. R. & Vine, S. G. 2003 The hierarchical formation of a stellar 

cluster. Mon. Not. R. Astron. Soc, 343, 413-418. 
Bonnell, I. A. & Bate, M. R. 2006 Star formation through gravitational collapse and 

competitive accretion. Mon. Not. R. Astron. Soc, 370, 488-494. 
Brodie, J. & Strader, J. 2006 Extragalactic globular clusters and galaxy formation, Annu. 

Rev. Astron. Astrophys., 44, 193-267. 
Bruzual A, G. 2010 Star clusters as simple stellar populations. Phil. Trans. R. Soc. A, 

this volume. 

Capuzzo Dolcetta, R., Di Matteo, P. & Miocchi, P. 2005 Formation and evolution of 

clumpy tidal tails around globular clusters. Astron. J., 129, 1906-1921. 
Carretta, E., et al. 2009a Na-0 anticorrelation and HB. VII. The chemical composition 

of first and second-generation stars in 15 globular clusters from GIRAFFE spectra. 

Astron. Astrophys., in press (arXiv:0909.2938 ). 
Carretta, E., Bragaglia, A., Gratton, R. & Lucatello, S. 20096 Na-O anticorrelation and 

HB. VIII. Proton-capture elements and metallicities in 17 globular clusters from UVES 

spectra. Astron. Astrophys., in press (arXiv:0909.2941 ). 
Cartwright, A. & Whitworth, A. P. 2004 The statistical analysis of star clusters. Mon. 

Not. R. Astron. Soc, 348, 589-598. 
Chernoff, D. F., Kochanek, C. S. & Shapiro, S. L. 1986 Tidal heating of globular clusters. 

Astrophys. J., 309, 183-209. 
Chernoff, D. F. & Shapiro, S. L. 1987 Globular cluster evolution in the Galaxy: a global 

view. Astrophys. J., 322, 113-122. 
Chernoff, D. F. & Djorgovski, S. 1989 An analysis of the distribution of globular clusters 

with postcollapse cores in the Galaxy. Astrophys. J., 339, 904-918. 
Chernoff, D. F. & Weinberg M. D. 1990 Evolution of globular clusters in the Galaxy. 

Astrophys. J., 351, 121-156. 
Clarke, C. 2010 The physics and modes of star cluster formation: simulations. Phil. Trans. 

R. Soc A, this volume. 



Article to appear in the Philosopical Transactions of the Royal Society A 



Star cluster dynamics 



15 



Cohn, H., Hut, P. & Wise, M. 1989 Gravothermal oscillations after core collapse in globular 

cluster evolution. Astrophys. J., 342, 814-822. 
Davis, D. S., Richer, H. B., Anderson, J., Brewer, J., Hurley, J., Kalirai, J. S., Rich, R. 

M. & Stetson, P. B. 2008a Deep Advanced Camera for Surveys imaging in the globular 

cluster NGC 6397: the binary fraction. Astron. J., 135, 2155-2162. 
Davis, D. S., Richer, H. B., King, I. R., Anderson, J., Coffey, J., Fahlman, G. G., Hurley, J. 

& Kalirai, J. S. 2008b On the radial distribution of white dwarfs in the globular cluster 

NGC 6397. Mon. Not. R. Astron. Soc, 383, L20-L24. 
Decressin, T., Meynet, G., Charbonnel, C, Prantzos, N. & Ekstrom, S. 2007 Fast rotating 

massive stars and the origin of the abundance patterns in galactic globular clusters. 

Astron. Astrophys., 464, 1029-1044. 
de Grijs, R. 2010 A revolution in star cluster research: setting the scene. Phil. Trans. R. 

Soc. A, this volume. 

de Grijs, R., Gilmore, G. F., Johnson, R. A. & Mackey, A. D. 2002 Mass segregation in 

young compact star clusters in the Large Magellanic Cloud. II. Mass functions. Mon. 

Not. R. Astron. Soc, 331, 245-258. 
Dehnen, W., Odenkirchen, M., Grebel, E. K. & Rix, H.-W. 2004 Modeling the disruption 

of the globular cluster Palomar 5 by Galactic tides. Astron. J., 127, 2753-2770. 
De Marchi, G., Paresce, F. & Pulone, L. 2007 Why haven't loose globular clusters collapsed 

yet? Astrophys. J., 656, L65-L68. 
Djorgovski, S. & King, I. 1986 A preliminary survey of collapsed cores in globular clusters. 

Astrophys. J., 305, L61-L65. 
Djorgovski, S. & Meylan, G. 1994 The Galactic globular cluster system. Astron. J., 108, 

1292-1311. 

D'Ercole, A., Vesperini, E., D'Antona, F., McMillan, S. L. W. & Recchi, S. 2008 Formation 

and dynamical evolution of multiple stellar generations in globular clusters. Mon. Not. 

R. Astron. Soc, 391, 825-843. 
Elmegreen, B. 2000 Star formation in a crossing time. Astrophys. J., 530, 277-281. 
Elson, R. A. W., Fall, S. M. & Freeman, K. C. 1987 The structure of young star clusters 

in the Large Magellanic Cloud. Astrophys. J., 323, 54-78. 
Ernst, A., Just, A. & Spurzem, R. 2009 On the dissolution of star clusters in the Galactic 

Centre. I. Circular orbits. Mon. Not. R. Astron. Soc, 399, 141-156. 
Fall, S.M. & Zhang, Q. 2001 Dynamical evolution of the mass function of globular star 

clusters. Astrophys. J., 561, 751-765. 
Fellhauer, M., Wilkinson, M. I. & Kroupa P. 2009 Merging time-scales of stellar subclumps 

in young star- forming regions. Mon. Not. R. Astron. Soc, 397, 954-962. 
Fischer, P., Pryor, C, Murray, S., Mateo, M. & Richtler, T. 1998 Mass segregation in 

young Large Magellanic Cloud clusters. I. NGC 2157. Astron. J., 115, 592-604. 
Fregeau, J. M. & Rasio F. A. 2007 Monte Carlo simulations of globular cluster evolution. 

IV. Direct integration of strong interactions. Astrophys. J., 658, 1047-1061. 
Fregeau, J. M., Richer, H. B., Rasio, F. A. & Hurley, J. R. 2009 The dynamical effects of 

white dwarf birth kicks in globular star clusters. Astrophys. J., 695, L20-L24. 
Fukushige, T. & Heggie, D. C. 1995 Pre-collapse evolution of galactic globular clusters. 

Mon. Not. R. Astron. Soc, 276, 206-218. 
Fukushige, T. & Heggie, D. C. 2000 The time-scale of escape from star clusters. Mon. Not. 

R. Astron. Soc, 318, 753-761. 
Gao, B., Goodman, J., Cohn, H. & Murphy, B. 1991 Fokker-Planck calculations of star 

clusters with primordial binaries. Astrophys. J., 370, 567-582. 
Giersz, M. & Heggie, D. C. 1994 Statistics of iV-body simulations. I. Equal masses before 

core collapse. Mon. Not. R. Astron. Soc, 268, 257-275. 



Article to appear in the Philosopical Transactions of the Royal Society A 



16 



E. Vesperini 



Giersz, M. & Spurzem, R. 2000 A stochastic Monte Carlo approach to model real star clus- 
ter evolution. II. Self-consistent models and primordial binaries. Mori. Not. R. Astron. 
Soc, 317, 581-606. 

Giersz, M. & Spurzem, R. 2003 A stochastic Monte Carlo approach to model real star 

cluster evolution. III. Direct integration of three- and four-body interactions. Mon. 

Not. R. Astron. Soc, 343, 781-795. 
Gieles, M. & Baumgardt, H. 2008 Lifetimes of tidally limited star clusters with different 

radii. Mon. Not. R. Astron. Soc, 389, L28-L32. 
Gill, M., Trenti, M., Miller, M. C, van der Marel, R., Hamilton, D. & Stiavelli, M. 2008 

Intermediate-mass black hole induced quenching of mass segregation in star clusters. 

Astrophys. J., 686, 303-309. 
Glebbeek, E., Gaburov, E., de Mink, S. E., Pols, O. R. & Portegies Zwart, S. F. 2009 The 

evolution of runaway stellar collision products. Astron. Astrophys., 497, 255-264. 
Gnedin, O. Y. & Ostriker, J. P. 1997 Destruction of the Galactic globular cluster system. 

Astrophys. J., 474, 223-255. 
Gnedin, O. Y., Lee, H. M. & Ostriker, J. P. 1999 Effects of tidal shocks on the evolution 

of globular clusters. Astrophys. J., 522, 935-949. 
Goodman, J. & Hut, P. 1989 Primordial binaries and globular cluster evolution. Nature, 

339, 40-42. 

Goodwin, S. P. 1998 Constraints on the initial conditions of globular clusters. Mon. Not. 

R. Astron. Soc, 294, 47-60. 
Goodwin, S. P. 2010 Binaries in star clusters and the origin of the field stellar population. 

Phil. Trans. R. Soc A, this volume. 
Goodwin, S. P. & Whitworth, A. P. 2004 The dynamical evolution of fractal star clusters: 

the survival of substructure. Astron. Astrophys., 413, 929-937. 
Gouliermis, D., Keller, S. C, Kontizas, M., Kontizas, E. & Bellas- Velidis, I. 2004 Mass 

segregation in young Magellanic Cloud star clusters: four clusters observed with HST. 

Astron. Astrophys., 416, 137-155. 
Gouliermis, D. A., de Grijs, R. & Xin, Y. 2009 Astrophys. J., 692, 1678-1689. 
Gratton, R., Sneden, C. & Carretta, E. 2004 Abundance variations within globular clus- 
ters. Annu. Rev. Astron. Astrophys., 42, 385-440. 
Grillmair, C. J. & Dionatos, O. 2006 A 22° tidal tail for Palomar 5. Astrophys. J., 641, 

L37-L39. 

Giirkan, M. A., Freitag, M. & Rasio, F. A. 2004 Formation of massive black holes in dense 
star clusters. I. Mass segregation and core collapse. Astrophys. J., 604, 632-652. 

Harris, W. E. 2001 Globular cluster systems. In Star Clusters (eds L. Labhardt & B. 
Binggeli), Saas-Fee Adv. Course 28, pp. 223-294, Berlin: Springer. 

Harris, W. E. 2010 Massive star clusters in galaxies. Phil. Trans. R. Soc A, this volume. 

Heggie, D. C. 1975 Binary evolution in stellar dynamics. Mon. Not. R. Astron. Soc, 173, 
729-787. 

Heggie, D. C. 1992 Ecology of globular clusters. Nature, 359, 772-773. 

Heggie, D. C. 2001 Mass loss from globular clusters. In Dynamics of Star Clusters and the 

Milky Way (eds S. Deiters, B. Fuchs, R. Spurzem, A. Just & R. Wielen), Astron. Soc. 

Pac. Conf. Ser., 228, 29-42. San Francisco: Astron. Soc. Pac. 
Heggie, D. C. & Hut, P. 1993 Binary-single-star scattering. IV. Analytic approximations 

and fitting formulae for cross sections and reaction rates. Astrophys. J. Suppl. Ser., 85, 

347-409. 

Heggie, D. C. & Aarseth, S. J. 1992 Dynamical effects of primordial binaries in star 
clusters. I. Equal masses. Mon. Not. R. Astron. Soc, 257, 513-536. 

Heggie, D. C. & Hut, P. 2003 The Gravitational Million-Body Problem, Cambridge: Cam- 
bridge University Press. 



Article to appear in the Philosopical Transactions of the Royal Society A 



Star cluster dynamics 



17 



Heggie, D. C, Trenti, M. & Hut, P. 2006 Star clusters with primordial binaries. I. Dy- 
namical evolution of isolated models. Mon. Not. R. Astron. Soc, 368, 677-689. 

Heggie, D. C, Hut, P., Mineshige, S., Makino, J. & Baumgardt, H. 2007 The core radius 
of a star cluster containing a massive black hole. Publ. Astron. Soc. Jpn, 59, L11-L14. 

Heggie, D. C. & Giersz, M. 2008 Monte Carlo simulations of star clusters. V. The globular 
cluster M4. Mon. Not. R. Astron. Soc, 389, 1858-1870. 

Heggie, D. C. & Giersz, M. 2009a Monte Carlo simulations of star clusters. VI. The 
globular cluster NGC 6397. Mon. Not. R. Astron. Soc, 395, 1173-1183. 

Heggie, D. C. & Giersz, M. 20096 1 Gyr in the life of the globular cluster NGC 6397. Mon. 
Not. R. Astron. Soc, 397, L46-L50. 

Henon, M. 1961 Sur revolution dynamique des amas globulaires. Ann. d'Astrophys., 24, 
369-419. 

Heyl, J. 2008 Constraining white-dwarf kicks in globular clusters. III. Cluster heating. 

Mon. Not. R. Astron. Soc, 390, 622-624. 
Heyl, J. & Penrice, M. 2009 Constraining white dwarf kicks in globular clusters. IV. 

Retarding core collapse. Mon. Not. R. Astron. Soc, 397, L79-L81. 
Hillenbrand, L. A. 1997 On the stellar population and star-forming history of the Orion 

Nebula Cluster. Astron. J., 113, 1733-1768. 
Hillenbrand, L. A. & Hartmann, L. E. 1998 A preliminary study of the Orion Nebula 

Cluster structure and dynamics. Astrophys. J., 492, 540-553. 
Hills, J. 1980 The effect of mass loss on the dynamical evolution of a stellar system. 

Analytic approximations. Astrophys. J., 235, 986-991. 
Hurley, J. 2007 Ratios of star cluster core and half-mass radii: a cautionary note on 

intermediate-mass black holes in star clusters. Mon. Not. R. Astron. Soc, 379, 93-99. 
Hurley, J. R. & Bekki, K. 2008 Self-consistent simulations of star cluster formation from 

gas clouds under the influence of galaxy-scale tidal fields. Mon. Not. R. Astron. Soc, 

389, L61-L65. 

Just, A., Berczik, P., Petrov, M. I. & Ernst, A. 2009 Quantitative analysis of clumps in 

the tidal tails of star clusters. Mon. Not. R. Astron. Soc, 392, 969-981. 
Kalirai, J. S. & Richer, H. B. 2010 Star clusters as laboratories for stellar and dynamical 

evolution. Phil. Trans. R. Soc A, this volume. 
Klessen, R. 2001 The formation of stellar clusters: mass spectra from turbulent molecular 

cloud fragmentation. Astrophys. J., 556, 837-846. 
Kumar, M. S. N. & Schmeja, S. 2007 The spatial distribution of substellar objects in IC 

348 and the Orion Trapezium cluster. Astron. Astrophys., 471, L33-L36. 
Kiipper, A. H. W., MacLeod, A. & Heggie, D. C. 2008 On the structure of tidal tails. 

Mon. Not. R. Astron. Soc, 387, 1248-1252. 
Koch, A., Grebel, E. K., Odenkirchen, M., Martfnez-Delgado, D. & Caldwell, J. A. R. 

2004 Mass segregation in the globular cluster Palomar 5 and its tidal tails. Astron. J., 

128, 2274-2287. 

Kroupa, P., Aarseth, S. & Hurley, J. 2001 The formation of a bound star cluster: from the 
Orion Nebula Cluster to the Pleiades. Mon. Not. R. Astron. Soc, 321, 699-712. 

Kruijssen, J. M. D. 2008 Explaining the mass-to-light ratios of globular clusters. Astron. 
Astrophys., 486, L21-L24. 

Kruijssen, J. M. D. & Mieske, S. 2009 Dissolution is the solution: on the reduced mass- 
to- light ratios of Galactic globular clusters. Astron. Astrophys., 500, 785-799. 

Krumholz, M. R., Klein, R. I. & McKee, C. F. 2005 The formation of stars by gravitational 
collapse rather than competitive accretion. Nature, 438, 332-334. 

Krumholz, M. R. & Bonnell, I. A. 2007 Models for the formation of massive stars. In 
Structure Formation in the Universe (ed. G. Chabrier), in press, Cambridge: Cambridge 
University Press l|arXiv:0712.0828|l 



Article to appear in the Philosopical Transactions of the Royal Society A 



18 



E. Vesperini 



Lada, C. J. 2010 The physics and modes of star cluster formation: observations. Phil. 

Trans. R. Soc. A, this volume. 
Lamers, H. J. G. L. M., Gieles, M., Bastian, N., Baumgardt, H., Kharchenko, N. V. & 

Portegies Zwart, S. 2005 An analytical description of the disruption of star clusters 

in tidal fields with an application to Galactic open clusters. Astron. Astrophys., 441, 

117-129. 

Lamers, H. J. G. L. M., Anders, P. & de Grijs, R. 2006 The photometric evolution of 
dissolving star clusters. I. First predictions. Astron. Astrophys., 452, 131-140. 

Larsen, S. S. 2004 The structure and environment of young stellar clusters in spiral galax- 
ies. Astron. Astrophys., 416, 537-553. 

Larsen, S. S. 2010 Young and intermediate-age massive star clusters. Phil. Trans. R. Soc. 
A, this volume. 

Leon, S., Meylan, G. & Combes, F. 2000 Tidal tails around 20 Galactic globular clusters. 
Observational evidence for gravitational disk/bulge shocking. Astron. Astrophys., 359, 
907-931. 

Lynden-Bell, D. & Wood, R. 1968 The gravo-thermal catastrophe in isothermal spheres 
and the onset of red-giant structure for stellar systems. Mon. Not. R. Astron. Soc, 138, 
495-525. 

Mackey, A. D., Wilkinson, M. L, Davies, M. B. & Gilmore, G. F. 2007 The effect of 

stellar-mass black holes on the structural evolution of massive star clusters. Mon. Not. 

R. Astron. Soc, 379, L40-L44. 
Mackey, A. D., Wilkinson, M. I., Davies, M. B. & Gilmore, G. F. 2008 Black holes and 

core expansion in massive star clusters. Mon. Not. R. Astron. Soc, 386, 65-95. 
McGlynn, T. A. 1984 Dissipationless collapse of galaxies and initial conditions. Astrophys. 

J., 281, 13-30. 

McLaughlin, D. E. & Fall, S. M. 2008 Shaping the globular cluster mass function by 

stellar-dynamical evaporation. Astrophys. J., 679, 1272-1287. 
McMillan, S. L. W. 2008 Gravitational dynamics of large stellar systems. Class. Quant. 

Grav., 25, 114007 (17pp). 
McMillan, S. L. W., Hut, P. & Makino, J. 1990 Star cluster evolution with primordial 

binaries. I. A comparative study. Astrophys. J., 362, 522-537. 
McMillan, S. L. W, Hut, P. & Makino, J. 1991 Star cluster evolution with primordial 

binaries. II. Detailed analysis. Astrophys. J., 372, 111-124. 
McMillan, S. L. W., Vesperini, E. & Portegies Zwart, S. F. 2007 A dynamical origin for 

early mass segregation in young star clusters. Astrophys. J., 655, L45-L49. 
Milone, A. P., et al. 2008 The ACS survey of Galactic globular clusters. III. The double 

subgiant branch of NGC 1851. Astrophys. J., 673, 241-250. 
Moeckel, N. & Bonnell, I. A. 2009 Does sub-cluster merging accelerate mass segregation 

in local star formation? Mon. Not. R. Astron. Soc, submitted (arXiv:0908.0253 ) 
Odenkirchen, M., Grebel, E. K., Dehnen, W., Rix, H.-W., Yanny, B., Newberg, H. J. 

Rockosi, C. M., Martmez-Delgado, D., Brinkmann, J. & Pier, J. R. 2003 The extended 

tails of Palomar 5: A 10° arc of globular cluster tidal debris. Astron. J., 126, 2385-2407. 
Odenkirchen, M., Grebel, E. K., Kayser, A. Rix, H.-W. & Dehnen, W. 2009 Kinematics 

of the tidal debris of the globular cluster Palomar 5. Astron. J., 137, 3378-3387. 
Parmentier, G., Goodwin, S. P., Kroupa, P. & Baumgardt, H. 2008 The shape of the 

initial cluster mass function: what it tells us about the local star formation efficiency. 

Astrophys. J., 678, 347-352. 
Pasquato, M., Trenti, M., De Marchi, G., Gill, M., Hamilton, D. P., Miller, M. C, Stiavelli, 

M. & van der Marel, R. P. 2009 Mass segregation in NGC 2298: limits on the presence 

of an intermediate mass black hole. Astrophys. J., 699, 1511-1517. 



Article to appear in the Philosopical Transactions of the Royal Society A 



Star cluster dynamics 



19 



Piotto, G. & Zoccali, M. 1999 HST luminosity functions of the globular clusters M10, 
M22, and M55. A comparison with other clusters. Astron. Astrophys., 345, 485-498. 

Piotto, G., Bedin, L. R., Anderson, J., King, I. R., Cassisi, S., Milone, A. P., Villanova, 
S., Pietrinferni, A. & Renzini, A. 2007 A triple main sequence in the globular cluster 
NGC 2808. Astrophys. J., 661, L53-L56. 

Portegies Zwart, S.F., Hut, P., Makino, J. & McMillan, S. L. W. 1998 On the dissolution 
of evolving star clusters. Astron. Astrophys., 337, 363-371. 

Portegies Zwart, S. F. & McMillan, S. L. W. 2002 The runaway growth of intermediate- 
mass black holes in dense star clusters. Astrophys. J., 576, 899-907. 

Portegies Zwart, S. F., Baumgardt, H., Hut, P., Makino, J. & McMillan, S. L. W. 2004 
Formation of massive black holes through runaway collisions in dense young star clusters 
Nature, 428, 724-726. 

Portegies Zwart, S .F., McMillan, S. L. W. & Makino, J. 2007 Star cluster ecology. VII. 

The evolution of young dense star clusters containing primordial binaries. Mori. Not. 

R. Astron. Soc, 374, 95-106. 
Renzini, A. 2008 Origin of multiple stellar populations in globular clusters and their helium 

enrichment. Mon. Not. R. Astron. Soc, 391, 354-362. 
Scally, A. & Clarke, C. 2002 Primordial substructure in the Orion Nebula Cluster. Mon. 

Not. R. Astron. Soc, 334, 156-166. 
Schmeja, S. & Klessen, R. S. 2006 Evolving structures of star-forming clusters. Astron. 

Astrophys., 449, 151-159. 
Sirianni, M., Nota, A., De Marchi, G., Leitherer, C. & Clampin, M. 2002 The low end of 

the initial mass function in young clusters. II. Evidence for primordial mass segregation 

in NGC 330 in the Small Magellanic Cloud. Astrophys. J., 579, 275-288. 
Stolte, A., Brandner, W., Brandl, B. & Zinnecker, H. 2006 The secrets of the nearest 

starburst cluster. II. The present-day mass function in NGC 3603. Astron. J., 132, 

253-270. 

Takahashi, K. & Portegies Zwart, S. F. 2000 The evolution of globular clusters in the 
Galaxy. Astrophys. J., 535, 759-775. 

Trenti, M. 2006 Dynamical evidence for intermediate-mass black holes in old globular 
clusters. Mon. Not. R. Astron. Soc, submitted (astro-ph/0612040). 

Trenti, M., Bertin, G. & van Albada, T. S. 2005 A family of models of partially relaxed 
stellar systems. II. Comparison with the products of collisionless collapse. Astron. As- 
trophys., 429, 161-172. 

Trenti, M., Heggie, D. C. & Hut, P. 2007a Star clusters with primordial binaries. II. 
Dynamical evolution of models in a tidal field. Mon. Not. R. Astron. Soc, 374, 344- 
356. 

Trenti, M., Ardi, E., Mineshige, S. & Hut, P. 20076 Star clusters with primordial binaries. 

III. Dynamical interaction between binaries and an intermediate-mass black hole. Mon. 

Not. R. Astron. Soc, 374, 857-866. 
Trenti, M., Vesperini, E. & Pasquato, M. 2009 Tidal disruption, global mass function and 

structural parameters evolution in star clusters, submitted, 
van Albada, T. S. 1982 Dissipationless galaxy formation and the R 1 ' 4 law. Mon. Not. R. 

Astron. Soc, 201, 939-955. 
van Loon, J. T. 2010 Chemical evolution of star clusters. Phil. Trans. R. Soc. A, this 

volume. 

Ventura, P., D'Antona, F., Mazzitelli, I. & Gratton, R. 2001 Predictions for self-pollution 

in globular cluster stars. Astrophys. J., 550, L65-L69. 
Vesperini, E. 1997 On the evolution of the Galactic globular cluster system. Mon. Not. R. 

Astron. Soc, 287, 915-928. 



Article to appear in the Philosopical Transactions of the Royal Society A 



20 



E. Vesperini 



Vesperini, E. 1998 Evolution ol the mass lunction of the Galactic globular cluster system. 

Mon. Not. R. Astron. Soc, 299, 1019-1039. 
Vesperini, E. 2000 Evolution of globular cluster systems in elliptical galaxies. I. Log- normal 

initial mass function. Mon. Not. R. Astron. Soc, 318, 841-856. 
Vesperini, E. 2001 Evolution of globular cluster systems in elliptical galaxies. II. Power-law 

initial mass function. Mon. Not. R. Astron. Soc, 322, 247-256. 
Vesperini, E. & Chernoff, D. F. 1994 The range of core size of postcollapse clusters sup- 
ported by primordial binaries. Astrophys. J., 431, 231-236. 
Vesperini, E. & Heggie, D. C. 1997 On the effects of dynamical evolution on the initial 

mass function of globular clusters. Mon. Not. R. Astron. Soc, 289, 898-920. 
Vesperini, E. & Zepf, S. E. 2003 Effects of the dissolution of low-concentration globular 

clusters on the evolution of globular cluster systems. Astrophys. J., 587, L97-L100. 
Vesperini, E., Zepf, S. E., Kundu, A. & Ashman, K. M. 2003 Modeling the dynamical 

evolution of the M87 globular cluster system. Astrophys. J., 593, 760-771. 
Vesperini, E., McMillan S. L. W. & Portegies Zwart S. 2009 Effects of primordial mass 

segregation on the dynamical evolution of star clusters. Astrophys. J., 698, 615-622. 
Weinberg, M. D. 1994 Adiabatic invariants in stellar dynamics. III. Application to globular 

cluster evolution. Astron. J., 108, 1414-1420. 
Zhang, Q. & Fall, S. M. 1999 The mass function of young star clusters in the Antennae 

galaxies. Astrophys. J., 527, L81-L84. 
Zepf, S. E. 2008 Observational constraints on the formation and evolution of globular 

cluster systems. In Dynamical Evolution of Dense Stellar Systems (eds E. Vesperini, 

M. Giersz & A. Sills), Proc. Int. Astron. Union Symp. 246, pp. 394-402, Cambridge: 

Cambridge University Press. 



Article to appear in the Philosopical Transactions of the Royal Society A 



